All sky CMB map from cosmic strings integrated Sachs- Wolfe effect 
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By actively distorting the cosmic microwave background (CMB) over our past light cone, cos- 
mic strings are unavoidable sources ol non-Gaussianity. Developing optimal estimators able to 
disambiguate a string signal from the primordial type ol non-Gaussianity requires calibration over 
synthetic lull sky CMB maps, which untill now had been numerically unachievable at the resolution 
ol modern experiments. In this paper, we provide the first high resolution lull sky CMB map ol 
the temperature anisotropies induced by a network ol cosmic strings since the recombination. The 
map has about 200 million sub-arcminute pixels in the healpix format which is the standard in use 
for CMB analyses (-/V a id e = 4096). This premiere required about 800,000 cpu hours; it has been 
generated by using a massively parallel ray tracing method piercing through thousands ol state ol 
art Nambu-Goto cosmic string numerical simulations which pave the comoving volume between the 
observer and the last scattering surface. We explicitly show how this map corrects previous results 
derived in the flat sky approximation, while remaining completely compatible at the smallest scales. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 

As all topological defects [l[ , cosmic strings incessantly 
generate gravitational perturbations all along the uni- 
verse history @-(3]- Their amplitude is directly given by 
the string energy density per unit length, GU <C 1 (in 
Planck unit), which is also the typical energy scale at 
which these objects are formed, eventually redshifted by 
warped extra-dimension for the so-called cosmic super- 
strings Although the theory of cosmological per- 
turbations can be applied to defects Q , predicting string 
induced anisotropies in the CMB is challenging. As op- 
posed to the perturbations of inflationary origin, which 
are generated once and for all in the early universe, active 
sources require the complete knowledge of their evolution 
at all times, from their formation till today [ltJ-Q- 

For these reasons, cosmological analyses often rely on 
analytical, or semi-analytical defect models [l5l - [25| which 
may not be accurate enough in view of the incoming 
flow of high precision CMB data, such as those from 
the Planck satellite and the other sub-orbital experi- 
ments [2t3 - [29j . The theoretical understanding of cos- 
mic string evolution in a Friedmann-Lemaitre- Robertson- 
Walker (FLRW) universe is still an active field of research 
which has led to the development of various theoretical 
models [3CH40J and numerical simulations, the latter hav- 
ing the advantage of incorporating all the defect dynam- 
ics [11-52]. However, within simulations, the dynami- 



* Electronic address: 
t Electronic address: 



christophc.ringcval@uclouvain.be 



bouchet@iap.fr 



cal range remains limited such that one has to extrapo- 
late numerical results over orders of magnitude by means 
of their scale invariant properties. The CMB tempera- 
ture and polarization angular power spectra have been 
derived for local strings only within the Abelian Higgs 
model [53U551 ] . As simulations do have only one free 
parameter that is GU, they provide a robust correspon- 
dence between string tension and CMB amplitude. From 
current data, Ref. [56j reports the two-sigma confidence 
limit GU < 4.2 x 10~ 7 , whereas semi-analytical methods 
find bounds ranging from to 10 -6 [57h60 |. 

Among other signatures, non-Gaussianities are un- 
avoidable consequences of the presence of cosmic strings 
(see Ref. [6l|, [62j for a review) . The determination of the 
power spectrum, i.e. the two-points function, from nu- 
merical simulations is not easy, and the situation is even 
worse for any higher n-point function. A way around 
this is to include photons inside a string simulation to 
produce a realization of the expected CMB tempera- 
ture anisotropies. In that respect, the resulting map 
contains all the statistical content, non-Gaussianities in- 
cluded. This method has originally been introduced for 
Nambu-Goto strings in Ref. [63[ and revived in Ref. [(54| 
to create a collection of statistically independent small 
angle CMB maps. As shown by Hindmarsh, Stebbins and 
Veeraraghavan, the small angle limit happens to be very 
convenient as the perturbed photon propagation equa- 
tions, namely the string induced integrated Sachs- Wolfe 
(ISW) effect [H, [6(| , reduces to a more tractable two- 
dimensional problem [g3, IHJ]. Those maps have been 
shown to be accurate as they correctly reproduce the 
small scale power spectrum of Abelian str ings f5oTl as 
well as the analytically expected one-point [69l - l7l| and 
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higher n-point functions [72H74j . Although fiat sky maps 
are adequate to devise new string-oriented searches in 
small scale CMB data [75l - [77l |. current searches for non- 
Gaussianities are mainly driven by the primordial type 
and based on full sky optimized estimators [781481] (see 
Ref. [HI for a review). 

In this paper, we generalize the method used in 
Ref. [64| and go beyond the flat sky approximation to 
generate a full sky CMB map induced by Nambu-Goto 
strings. The most unequivocal signature from strings are 
the temperature discontinuities they induce, which are 
naturally most striking at small scales; our efforts thus 
have been to achive a high resolution over the complete 
sky. In a hierarchical equal area isolatitude pixelisation 
of the sky, we have been able to maintain an angular res- 
olution of 8 les = 0.85', i.e. using the publicly available 
HEALPix code [Hj], our map has A side = 4096, i.e. 2 x 10 s 
pixels. As detailed in the following, our method includes 
all string effects from the last scattering surface till today, 
but does not include the Doppler contributions induced 
by the strings into the plasma prior to recombination. As 
a result, our map represents the ISW contribution from 
strings, which is dominant at small scales but underesti- 
mates the signal on intermediate length scales. Including 
the Doppler effects requires the addition of matter in the 
simulations, an approach which has been implemented in 
Refs. |84l85j | and recently used to generate a full sky map 
in Ref. [86|. As discussed in that reference, the comput- 
ing resources required to include matter severely limit the 
achievable resolution to 14' (N sidc = 256, with 800 000 
pixels). In that respect, our map is complementary to 
the one of Ref. [86} while extending the domain of appli- 
cability of existing small angle maps 64]. In particular, 
we recover the turnover in the spectrum observed around 
I ~ 200 in Ref. j86f] , which was cut by the small field of 
view of the flat sky maps. 

The paper is organized as follows. In Sec. [HI we briefly 
recap the characteristics of our numerical simulations 
of cosmic string networks and describe the ray tracing 
method used to compute the full sky map. The map it- 
self is presented in Sec. IIIII and compared with the small 
angle maps in the applicable limit. We also compute the 
angular power spectrum and conclude in the last section. 



II. METHOD 

A. All sky string ISW 

Denoting by X^(r, a) the string embedding functions, 
in the transverse temporal gauge 1 , up to a dipole term, 
the integrated Sachs- Wolfe contribution sourced by the 



Nambu-Goto stress tensor reads [68[ 
9(n) = -4GU f u ■ 



X 



: da, 



(1) 



where O(n) stands for the relative photon temperature 
shifts (8 = AT/T) in the observation direction n. The 
integral is performed over all string position vectors X — 
{X 1 } intercepting our past line cone, the invariant string 

length element being d^ = eda with e 2 = X /(l — X ). 
The string dynamical effects are encoded in 



X 



(h-X)-X 
1 + n-X ' 



(2) 



X ■ X = 0, X° = r] = t, r) is the conformal time. 



where the "acute accent" and the "dot" denote differ- 
entiation with respect to the world sheet coordinates a 
and r, respectively. The small angle and flat sky approx- 
imations consists in taking the limit Xn — > X which 
assumes that the observation direction matches with the 
strin g lo cation and that all photon trajectories are par- 
allel [H [H, i3 ■ Let us notice that Eq. (fTJ cannot be 
applied to a straight static string, but such a situation 
never occurs for the realistic string configurations studied 
in the following [68|, [87j . 

In the following, we use the all sky expression of Eq. (fTJ) 
to compute the induced temperature anisotropics 0(n) 
in each of the wanted 2 x 10 8 pixelized directions. Equa- 
tion JT]) shows that, in each direction, one has to sum 
up the contribution of all string segments d^ intercepting 
our past light cone since the last scattering surface and 
determine, for each of them, X and X. Although one 
string lying behind the observer does not contribute more 
than a few percent to the overall signal, it is impossible to 
artificially cut it without adding spurious discontinuities 
in the map, which would dangerously mimic real string 
patterns. As discussed below, we have filled the comov- 
ing volume between today and the last scattering sur- 
face with a few thousand Nambu-Goto numerical simula- 
tions in Friedmann-Lemaitre- Robertson- Walker (FLRW) 
space-time. Typically, hundreds of millions of string seg- 
ments intercept on our past light cone, each of them has 
to be included in Eq. ([IJ to get the overall signal for one 
pixel. One can immediately understand the computing 
challenge to obtain a full sky map as the total number of 
expected iterations roughly sums up to 10 16 . 



B. Nambu-Goto string simulations 

In order to get a realistic string configuration between 
the last scattering surface and today, we have followed 
Refs. [63|, [64j and stacked FLRW string simulations us- 
ing an improved version of the Bennett and Bouchet 
Nambu-Goto cosmic string code [43|, H^]. The runs have 
the same characteristics as those used in Ref. [iij and, 
in particular, we include only the loops having a size 
larger than a time-dependent cut-off. The reason being 
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that loops smaller than this cut-off have a distribution 
which is known to be contaminated by relaxation effects 
from the numerical initial conditions. As explained in 
Ref. ffrij , the cutoff is dynamically chosen by monitoring 
the time evolution of the energy density distribution as- 
sociated with loops of different sizes. This ensures that 
we include only loops having an energy density evolving 
as in scaling, i.e. in 1/t 2 . One may be worried about 
the deficit due to the missing loops artificially removed 
by the cut-off. An upper bound of the systematic er- 
rors that could be induced by these relaxation effects can 
be found in that same reference (see Sec. II. D). It does 
not exceed 10% on average and concerns only very small 
scales. In the following, we recap some of the relevant 
physical properties underlying our numerical simulations 
(see Sec. II. B in Ref. [64J for more details). Each simula- 
tion allows one to trace the time evolution of a network 
of cosmic strings in scaling over a cubic comoving volume 
of typical size 



30- 



V^rad/^matVl + z i 
W^matVl + z i 



Gpc. 



(3) 



In this expression, z\ is the redshift at which the sim- 
ulation is started, h is the reduced Hubble parameter 
today and f2 ma t) ^rad are the density parameters of mat- 
ter and radiation today. Starting at the last scattering 
surface, zj = 1089, gives a simulation comoving box of 
L s - lm ~ 1.7 Gpc 2 , which corresponds to an angular size 
of 7.2°. At the same time strings are evolved, we propa- 
gate photons along the three spatial directions and record 
all X and X for each string segment intercepting those 
light cones. Depending on the simulation realization, and 
its location, this corresponds to typically 10 4 -10 5 pro- 
jected string segments. One run is limited in time, as 
we use periodic boundary conditions, and ends after a 
30-fold increase in the expansion factor, i.e. at a red- 
shift z e = 36. As a result, covering the whole sky re- 
quires stacking side-by-side many different simulations, 
all starting at z\ = 1089 and ending at z c — 36. The 
missing redshift range can be dealt exactly in the same 
manner by stacking another set of runs which start at 
zi = 36 and end at z c = 0.2. From Eq. (|3J|, we see that 
the low-redshift simulations have a size of L s im — 13 Gpc 
such that only a few of them are required to cover the 
whole comoving volume. Notice that the use of different 
simulations to fill the comoving space does not induce 
visible artifacts in the final map. The signal being only 
sourced by the subset of string segments intercepting the 
past light cone, the probability of seeing an edge is al- 
most vanishing. Finally, as in Ref. 64j . we have skipped 
the last interval from z = 0.2toz = 0as almost no string 
intercepts our past line cone in that range. 



In the next section, we describe in more detail how we 
cut and stack the cosmic string simulations. 



C. Stacking simulations with healpix cones 

The hierarchical equal area iso-latitude pixelisation 
(healpix) of the two-sphere is a efficient method to pix- 
elize the sky which is well-suited to and commonly used 
for CMB analysis (83|. Our stacking method relies on 
the voxelization of the three-dimensional ball using cones 
subtending healpix pixels on its boundaries, i.e. for the 
two redshifts at which we start and stop the numeri- 
cal simulations. For the high-redshift contributions, the 
healpix cones fill the whole spherical volume in between 
the last scattering surface at z\ — 1089 and the two- 
sphere at z c = 36 (see Fig. [T]). The rest of the volume, 
associated with the low redshift simulations, can be vox- 
elized exactly in the same way but starting on the two- 
sphere at z\ = 36 and ending at z c = 0.2. The actual 
string dynamics simulations are evolved in cubic comov- 
ing boxes in which we take only strings living inside a 
healpix cone. In order to minimize the number of simu- 
lations required, the problem is now reduced to find the 
largest healpix cone fitting inside a cubic comoving box 
for all redshifts of interest. Moreover, one has to ensure 
that the photons intercepting the strings travel towards 
the observer. As seen in Fig.Q] both requirements can be 
implemented by adequately rotating the simulation box 
such that photons face the observer line of sight n, and 
the farthest healpix pixel fits inside the farthest squared 
face of the simulation box. As we propagate three photon 
waves in each simulation, one can use the same simula- 
tion rotated three times. Keeping only the strings living 
inside the healpix cones, one needs to patch them up 
till the whole comoving volume is filled. For this pur- 
pose, we have used the algorithms implemented within 
the HEALPix library [sH ], 

For the high-redshift contribution, the above require- 
ments are satisfied with a healpix voxelization having 
N*? de = 16, i.e. 3072 truncated cones from z = 1089 
to z = 36, therefore calling for 3072 cosmic string runs. 
As we can use one run three times, we have only per- 
formed 1024 independent simulations as described in 
Sec. Ill Bl fand Ref. [64]]). Concerning the low-redshift con- 
tribution, the simulation volume L s im being much larger, 
Nlf de = 1 or 12 simulations are enough to cover the vol- 
ume. Let us emphasize, at this stage, that the above- 
mentioned HEALPix resolutions only concern the simula- 
tion stacking method. In the next section, we discuss 
how we discretize the CMB sky using another, but more 
conventional, healpix scheme. 



2 In the runs, Z/ S im is computed exactly within the ACDM model 
whereas Eq. (|3]l is an analytical approximation assuming no cos- 
mological constant and n rat j -C Slmat- 



D. Healpix sky 

Assuming we have stacked the cosmic string simula- 
tions as previously explained, we have at our disposal a 
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FIG. 1: Voxelization of the comoving space between the two 
redshifts at which cosmic string simulations start and stop, 
here z\ = 1089 (outer sphere) and z c — 36 (inner sphere). All 
strings living in each healpix cone are kept and stacked to fill 
the whole comoving volume. Our high redshift voxelization 
scheme uses iV2| e = 16, i.e. requires 3072 string simulations 
while the low redshift one, starting at Zi — 36 and ending 
now, has A^f de = 1, i.e. 12 simulations. 



realization of all X and X lying on our past light cone 
since the last scattering surface. From Eqs. ([TJ and ©, 
the remaining step is to actually perform the integral for 
each desired value of the observer direction n. Choosing 
the values of n has been made by using another healpix 
pixelization scheme, this time on the simulated sky. A 
typical Planck-like CMB experiment requiring an angu- 
lar resolution of 5', we have chosen an angular resolution 
of 0.9' to reduce sufficiently the small scale inaccuracies 
which are present at the subpixel scale. The correspond- 
ing healpix resolution is N s n c = 4096, i.e. a sky map 
having 2 x 10 8 pixels. 

In the next section, after having briefly exposed how 
the above computing challenges have been solved, we 
present the simulated CMB map. 

III. STRING SKY MAP 

The method we have exposed in the previous section 
has the advantage to be completely factorizable into two 
independent numerical problems. The first is to perform 
string simulations and record only those events intercept- 
ing the light cones. The second is to use these events to 
actually perform the integration of Eq. ([1]), which gives 
the final signal 8(n). 

A. Computing resources 

Performing the one thousand cosmic string runs has re- 
quired around 300,000 cpu- hours on current x86-64 pro- 
cessors. The needed memory and disk space resources 
remaining reasonable, the computations have used local 



computing resources provided by the Planck-HFI com- 
puting centre at the "Institut d'Astrophysique de Paris" 
and the CP3 cosmo cluster at Louvain University. By 
the end of the runs, the total number of string segments 
recorded on the light cone account for typically 1 Tb of 
storage data. 

From the light cone data, performing 2 x 10 8 line of 
sight integrals over 10 8 string segments [see Eq. (JXJ) ] is 
a serious numerical problem. This part of the code has 
therefore been parallelized at three levels using the sym- 
metries of the problem. 

First, we have used a distributed memory paralleliza- 
tion as implemented in the message passing interface 
(MPI) to split the string contributions into sub-blocks. 
From our previous discussing, a natural implementation 
is an MPI-parallelization over the healpix cones which 
therefore allows different machines to deal with a subset 
of cones. More intuitively, it means that the final map 
is obtained by adding up full sky maps, each one being 
sourced by the strings contained in a few healpix cones 
only. Second, for each of those processes, the compu- 
tation of the 2 x 10 8 pixels has been parallelized using 
the shared memory OpenMP directives. In other words, 
pixels can be simultaneously computed using all of the 
available processors inside a single machine. Finally, for 
each of the above OpenMP threads, we have vectorized the 
most inner loop, i.e. the discrete version of Eq. ([I]). Do- 
ing this allows to use simultaneously multiple registers of 
each processing unit to add up a few string segments at 
once. 

The computing resources have been provided by the 
National Energy Research Scientific Computing Center 
(NERSC) at the Lawrence Berkeley National Labora- 
tory 3 . The run has been performed on the hopper super- 
computer using 512 MPI nodes at the first level of paral- 
lelization. On each node, we have deployed the OpenMP 
parallelization over the 24 cores available. Each node 
is made of two twelve-cores "AMD MagnyCours" cpus 
which supports only a limited amount of vectorization. 
However, a vectorization over 16 string segments signifi- 
cantly improves cache memory latency and gives a final 
speed-up of two compared to a pure scalar processing. 
All in all, the whole computation required 12, 000 cores 
and has been completed after 500, 000 cpu-hours. 



B. Results 

In Fig. [21 we have plotted the final map together with 
the two contributions coming from high and low redshift. 
As expected, most of the strings are at high redshift, 
the low redshift contribution showing only a few strings 
crossing the sky. For this reason, we have not included 
strings with z < 0.2, as almost no strings are present. 
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String ISW 



100.0 DT/T/GU 



High Redshift Contribution 



100.0 DT/T/GU 



Low Redshift Contribution 




100.0 DT/T/GU 



FIG. 2: All sky CMB map of string induced ISW effect (upper panel) in a Mollweide projection. The two lower maps represent 
respectively the high redshift contribution (from last scattering to z = 36) and the low redshift contribution (from z = 36 to 
now). The color scale indicates the range of Q/GU fluctuations which, once multiplied by the value of GU and Tcmb, would 
thus be in the tens of (jlK range for GU ~ 10 -6 . 
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At first glance these maps may look Gaussian, which is 
because the string patterns essentially show up at small 
scales while they are averaged on the largest angular 
scales 86] . In Fig. [3J we have represented a zoom over a 
7.2° region in which one recovers the same string discon- 
tinuities as p reviously derived within the flat sky approx- 
imation |64j [. In order to make the comparison sharper, 
both the gnomic projection of our spherical patch and the 
flat sky map coming from the same string simulation are 
represented in Fig. [3J Up to some spherical distortions 
off-center, both maps predict the same structures. 

Finally, in Fig. 01 we have plotted the angular power 
spectrum obtained from the full sky map. In this figure, 
the dashed curve represents the mean value previously 
derived using flat sky patches in Ref . [64] . We observe a 
very small loss of power compared to the flat sky result, 
which may not be significant as it does not exceed the 
one-sigma fluctuations expected between various string 
realizations. However, as Fig. [3J suggests, the flat sky 
maps are a bit sharper than the spherical ones due to 
some missing spherical effects and, as such, they may 
slightly overestimate the signal. This is not surprising 
as the flat sky maps are derived using the small angle 
approximated version of Eq. ((TJ) where all scalar products 
are postulated to be either parallel or null [67J. As our 
full sky map does not approximate anything, it should 
contain slightly less power. 

Let us also notice that there is some extra power for I > 
5000 associated with the full sky power spectrum. This is 
a spurious aliasing effect coming from the slow decrease 
of the power spectrum at small scales combined with our 
ray tracing method. Each pixel of the map represents the 
real signal, but only in the centroid direction fi. This has 
the effect of including string structures smaller than the 
pixel angular resolution thereby aliasing the final map. 

In Fig.SJ the power spectrum turns-over for I < 200 as 
those modes correspond to wavelengths larger than the 
correlation length at recombination. This property was 
also observed in the map derived in Ref. [86[ but barely 
visible in the flat sky maps of Ref. due to their small 
field of view. We see, a posteriori, that this effect was 
however indeed present as the dash curve of Fig.|4]flattens 
at the same location than the full sky spectrum. 



IV. CONCLUSION 

The main result of this article is the map displayed in 
Fig. [2] which provides a realization of the all sky CMB 
temperature anisotropics induced by a network of cos- 
mic strings since the last scattering surface. The chal- 
lenges underlying this map come from both the require- 
ment of covering the whole sky and an unprecedented 
angular resolution of 0.85 minute of arc, associated with 
a HEALPix resolution of = 4096. However our map 
includes only the ISW string contribution. This is the 
dominant signal at small scales, but it misses the Doppler 
effects around the intermediate multipoles. As those ef- 
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70.0 DT/T/GU 



Flat sky 
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FIG. 3: Comparison between full sky and flat sky / small an- 
gle approximation. The upper panel is a gnomic projection 
of a 7.2° patch cut from the high redshift full sky map at 
the north pole. The lower panel is from Ref. [6J] and repre- 
sents the high redshift flat sky calculation coming from the 
same string simulation. Both maps exhibit globally the same 
anisotropy patterns and amplitude. More precisely, there is 
essentially no difference in the center, whereas geometrical 
and amplitude distortions become increasingly apparent to- 
wards the edges. This is expected as spherical effects are not 
included in the flat sky approximation. 
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fects have been computed in Ref. [86|, at the expense of 
having a poor resolution, it will be interesting to inves- 
tigate whether both results can be combined to obtain a 
fully accurate representation of the stringy sky over the 
full range of observable scales. 
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FIG. 4: Angular power spectrum extracted from the string 
map of Fig. [5] (solid red curve) compared to the one estimated 
in Ref. [6J], using the small angle and flat sky approximation 
(dashed black curve). The shaded region represents the one- 
sigma fluctuations over various string realizations (flat sky). 
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